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| „ , Abstract 

The subject matter of this paper concerns anisotropic diffusion equations: we con- 
£C) , sider heat equations whose diffusion matrix have disparate eigenvalues. We determine 
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first and second order approximations, we study the well-posedness of them and estab- 



p I ■ lish convergence results. The analysis relies on averaging techniques, which have been 

used previously for studying transport equations whose advection fields have disparate 



components. 
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1 Introduction 

Many real life applications lead to highly anisotropic diffusion equations: flows in porous me- 
dia, quasi-neutral plasmas, microscopic transport in magnetized plasmas [6], plasma thrusters, 
image processing [10], [12], thermal properties of crystals [9j. In this paper we investigate the 



behavior of the solutions for heat equations whose diffusion becomes very high along some 
direction. We consider the problem 

d t u £ - dw y (D(y)V y u £ ) - ^dw y (b(y) %)V/) = 0, (t, y)eR+x R m (1) 

u s (0,y)=uf n (y), y£R m (2) 

where D(y) £ A4 m (M) and b(y) € M. m are smooth given matrix field and vector field on M. m , 
respectively. For any two vectors £, rj, the notation ^ i] stands for the matrix whose entry 
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is and for any two matrix A, B the notation A : B stands for trace(*^4-B) = AijBij 
(using Einstein summation convention). We assume that at any y G M m the matrix D{y) is 
symmetric and D(y) + b(y) <S> b(y) is positive definite 

= D(y), 3d>0 such that D(y)£ ■ £ + (b(y) ■ £) 2 > d |£| 2 , (GR m , y G K m . (3) 

The vector field to which the anisotropy is aligned, is supposed divergence free i.e., 

divyb = 0. We intend to analyse the behavior of <(TJ), (J2J) for small e, let us say < e < 1. 
In that cases D(y) + \b(y) £>$ b(y) remains positive definite and if (uf n ) £ remain in a bounded 
set of L 2 (M m ), then (u £ ) £ remain in a bounded set of L°°(IR + ; L 2 (M. m )) since, for any t G 
we have 

\f (u £ (t,y)fdy + d f! \V y u e {s,y)\ 2 dyds<\ [ (u £ (t, y)) 2 dy 
z Jr™ Jo JR m z JM. m 

+ [[ \D(y) + -b(y)®b(y)\ : V y u £ (s,y) ®V y u £ (s,y) dyds 
JO JR m I £ ) 

z Jm m 

In particular, when e \ 0, (u £ ) e converges, at least weakly * in L°°(R + ; L 2 (M m )) towards 
some limit u G L 2 (M m )). Notice that the explicit methods are not well adapted for 

the numerical approximation of ([T]), d2D when e \ 0, since the CFL condition leads to severe 
time step constraints like 

d At 1 

IJKyf ~ 2 

where At is the time step and Ay is the grid spacing. In such cases implicit methods are 
desirable [I], [II]. 

Rather than solving ([1]), (J2J) f° r small e > 0, we concentrate on the limit model satisfied by 
the limit solution u = lim^o u £ . We will see that the limit model is still a parabolic problem, 
decreasing the L 2 (M m ) norm and satisfying the maximum principle. At least formally, the 
limit solution u is the dominant term of the expansion 

u £ = u + eu 1 +e 2 u 2 + ... (4) 

Plugging the Ansatz into ([I]) leads to 

div y {b ® VSJyu) = 0, (t, y) G K+ x R m (5) 

d t u-divy(DVyu) -diYyib^bVyU 1 ) = 0, (t,y) el + xl ra (6) 
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Clearly, the constraint ([5]) says that at any time t G M + , b ■ V y u = 0, or equivalently u(t, 
remains constant along the flow of b, see (fT5j) 



u(t,Y(s;y)) = u(t,y), seR, y G 



The closure for u comes by eliminating u 1 in (0) , combined with the fact that (0) holds true at 
any time t G R+. The symmetry of the operator div y (b®bV y ) implies that dtu — div y (DX7 y u) 
belongs to (ker(6 • V y ))^ and therefore we obtain the weak formulation 



d 
dt 



[ u{t,y)v(y)dy+ I DV y u(t,y)-Vy<p(y)dy = 0, ^ G H 1 (R m ) n ker(6 • V y ) . (7) 



The above formulation is not satisfactory, since the choice of test functions is constrained 
by ©; (0) is useless for numerical simulation. A more convenient situation is to reduce (J7J) 
to another problem, by removing the constraint ©. The method we employ here is related 
to the averaging technique which has been used to handle transport equations with diparate 
advection fields [2], [3], @], [5] 

d t u £ + a(t, y) ■ V y u £ + -b{y) ■ V y u £ = 0, (t, y) G M+ x R m (8) 

u e (0,y)=uf Q (y), yeR m . (9) 

Using the same Ansatz (jU) we obtain as before that b ■ V y u(t, •) = 0, t G M + and the closure 
for u writes 

Pr ojker(6-V w ){^ u + °- V y u} = (10) 

or equivalently 

— / u(t,y)(p(y) dy - f u(t, y) a ■ V y p dy = (11) 

for any smooth function satisfying the constraint b-V y ^p = 0. The method relies on averaging 
since the projection on ker (b- V^) coincides with the average along the flow of b, cf. Proposition 
13.11 As u satisfies the constraint b ■ V y u = 0, it is easily seen that Proj ker ( b . v y )dtu = dtu. 
A simple case to start with is when the transport operator a ■ V y and b ■ V y commute i.e., 
[b ■ Vy, a ■ Vy] = 0. In this case a ■ V y leaves invariant the subspace of the constraints, implying 
that Proj ker ( b . Va ){a • S/ y u} = a ■ V y u. Therefore (fTU|) reduces to a transport equation and 
it is easily seen that this equation propagates the constraint, which allows us to remove it. 
Things happen similarly when the transport operators a ■ V y , b ■ V y do not commute, but the 
transport operator of the limit model may change. In [3] we prove that there is a transport 
operator A ■ Vy, commuting with b ■ V y , such that for any u G ker (b ■ V y ) we have 

Proj ke r(&-V H ){ a ■ V a «} = A ■ V y U. 
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Once we have determined the field A, (|10p can be replaced by dtu + A ■ V y u = 0, which 
propagates the constraint b ■ V y u{t) = as well. 

Comming back to the formulation ((7]), we are looking for a matrix field D(y) such that 
divy(D'Vy) commutes with b ■ V y and 

P™ikcr(b-\7 y ){divy{D(y)Vyu)} = div y (D(y)V y u), u G ker(6 • V„). 

We will see that, under suitable hypotheses, it is possible to find such a matrix field D, and 
therefore ([7]) reduces to the parabolic model 

d t u - divy(D(y)X7 y u) = 0, (i, y) G R+ x R m . (12) 

The matrix field D will appear as the orthogonal projection of the matrix field D (with 
respect to some scalar product to be determined) on the subspace of matrix fields A satisfying 
[b • Vy, divylAVy)} = 0. The field D inherits the properties of D, like symmetry, positivity, 
etc. 

Our paper is organized as follows. The main results are presented in Section [2 Section 
El is devoted to the interplay between the average operator and first and second order linear 
differential operators. In particular we justify the existence of the averaged matrix field D 
associated to any field D of symmetric, positive matrix. The first order approximation is 
justified in Section H] and the second order approximation is discussed in Section [SJ Several 
technical proofs are gathered in Appendix [XI 

2 Presentation of the models and main results 

We assume that the vector field b : R m — > M m is smooth and divergence free 

b e W^(R m ), div,6 = (13) 

with linear growth 

3C>0 such that \b(y)\ < C(l + \y\), y G M. m . (14) 

We denote by Y(s;y) the characteristic flow associated to b 

^- = b(Y(s;y)), Y(s;0)=y, sCl, y G M m . (15) 

Under the above hypotheses, this flow has the regularity Y G W^q^°(M x R m ) and is measure 
preserving. 
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We concentrate on matrix fields A(y) G Ll oc (R m ) such that [b(y) ■ V y , dw y (A(y)Vy)] = 0, 
let us say in V(W n ). We check that the commutator between b ■ V y and dYVy{AV y ) writes 
cf. Proposition 13.71 

[b(y) ■ V y ,divy(A(y)S/ y )] = div„([6, A]V V ) in V'(R m ) 

where the bracket between b and A is given by 

[b, A] := (b ■ V y )A - d y bA{y) - A(y) %b, y € R m . 

Several characterizations for the solutions of [b, A] = in T> '(R m ) are indicated in the Propo- 
sitions 13.81 13.91 among which 

A(Y(s;y)) = dyY(s;y)A(y) t dyY(s;y), s£l, y € M. m . (16) 

We assume that there is a matrix field P(y) such that 

t P = P, P(y)£-Z>0, £eK m 5 yeR m , P" 1 , P £ Xf oc (M m ), [6, P] = in 2? '(M m ). (17) 

We introduce the set 

H Q = {A = A{y) : / Q{y)A{y) : A{y)Q(y) dy < +00} 

where Q = P _1 , and the scalar product 

(A,B) Q = / QA:BQdy, A,BeH Q . 

The equality (fT6j) suggests to introduce the family of applications G(s) : Hq — > Pq, s£R, 
= (<9 y Y) _1 (,s; •)^4(Y'(s; •)) t (d y Y)~ 1 (s; •) which is a C°-group of unitary operators 
on Hq cf. Proposition 13.121 This allows us to introduce L, the infinitesimal generator of 
(Gr(s)) S £R. The operator L is skew-adjoint on Hq and its kernel coincides with {A € Pq C 
L ioc( Rm ) '■ l b > A ] = in V'(R m )} cf. Proposition [3T3J The averaged matrix field denoted 
(D)q, associated to any D G Pq appears as the long time limit of the solution of 

d t A - L(L(A)) = 0, teR+ (18) 

A(0) = P. (19) 
The notation (•) stands for the orthogonal projection (in L 2 (M m )) on ker(6 • Vy). 
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Theorem 2.1 Assume that (fig]). (fTijl . (fTT]) reo/d frue. Then for any D G ifQnL°°(]R m ) £/ie 
solution of (fT8|) . (fT9j) converges weakly in Hq as t — ► +00 towards the orthogonal projection 
of D on ker L 

lim = <D)o wea£% in iT Q , <D) q := Proj kcri £>. 



If D is symmetric and positive, then so is the limit {D)q = limt_>+oo A(t), and satisfi 



es 



L((D) Q ) = 0, V y u- (D) Q V y v = {S/yU- DVyv) , u,v E H 1 (R m ) n ker (6 • V y ) (20) 



V y u- (D) Q V y (b- V^V)/ =0, u G il (R m ) n ker (6 • V y ), ip G C,r(lR m ). (21) 
The first order approximation (for initial data not necessarily well prepared) is justified by 



Theorem 2.2 Assume that (|13p . flHJ) , (|17p . (|5ip /10/d frae and that D is afield of symmetric 
positive matrix, which belongs to Hq. Consider a family of initial conditions (uf n ) e C L 2 (R m ) 
such that ({uf n )) £ converges weakly in L 2 (R m ), as e \ 0, towards souic function lijn- 
denote by u £ the solution of ([I]), ([2]) and n free solution of 

d t u-divy{(D) Q V y u) = 0, teR + , t/£l ra (22) 

u(0,y)=«to(lO, ^ r (23) 
where (D)q is associated to D, cf. Theorem \2.1\ Then we have the convergences 

limu e = u weakly* in L°°(R + ; L 2 (R m )) 

lim V y u £ = V y u weakly in L 2 (R + ; L 2 (R m )). 

The derivation of the second order approximation is more complicated and requires the 
computation of some other matrix fields. For simplicity, we content ourselves to formal 
results. The crucial point is to introduce the decomposition given by 

Theorem 2.3 Assume that (|14|) . (|17|) . (|51|) hold true and that L has closed range. 

Then, for any field of symmetric matrix D G Hq, there is a unique field of symmetric matrix 
F G dom(L 2 ) n (kerL)^ such that 

-diVy(DVy) = -diVy((D) Q Vy)+dWy(L 2 (F)Vy) 
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that is 

/ DVyU ■ V y v dy — {D)q V y u ■ V y v dy 

jR m ' JR m 

= / L{F)V y u-Vy{b -V y v) dy + / L(F)V y (b ■ W y u) ■ V y v dy 

jR m JR" 1 

= - FV y (b ■ V y (b ■ Vyu)) ■ Vyv dy-2 FV y (b ■ V y u) ■ V y (b ■ V y v) dy 

jR m JR m 

- I FVyU-Vy{b-Vy(b-VyV))dy 

jR m 

for any u,v G C^(R m ). 

After some computations we obtain, at least formally, the following model, replacing the 
hypothesis (|17p by the stronger one: there is a matrix field R(y) such that 

det R(y) ^ 0, y G R m , Q = t RR and P = Q^ 1 G Lf oc (R m ), b ■ W y R + Rd y b = in £> '(R m ). 

(24) 

Theorem 2.4 Assume that ([13]), CUD, ([28]), ([EE]), flM]) hold true and that D is afield of sym- 
metric positive matrix which belongs to HQnL° (R m ). Consider a family of initial conditions 
(n? n ) e C L 2 (R m ) such that ( " m ) £> o converges weakly in L 2 (R m ), as e \ 0, towards a 
function v- m , for some function u- in G ker(6- Vy). Then, a second order approximation for ([1]) 
is provided by 

d t u £ -div y ({D) Q V y u £ )+e[dwy{{D) Q Vy),div y (FV y )}u £ -eS(u £ ) = 0, (t,y) G R+ x R m 

(25) 

u £ (0, y) = u- m (y) + e(v in (y) + w in (y)), w in = dw y (FVyU in ), y G R m (26) 

for some fourth order linear differential operator S, see Proposition \5.!A and the matrix field 
F given by Theorem\2. 



3 The average operator 

We assume that the vector field b : R m —> R m satisfies (|13p . (j!4|) . We consider the linear 
operator u — > b ■ W y u = div y (u6) in L 2 (R m ), whose domain is defined by 

dom(6 • V y ) = {u G L 2 (R m ) : div y (ub) G L 2 (R m )}. 

It is well known that 

ker(6 • V„) = {u G L 2 (R m ) : u(Y(s; •)) = u(-), s G R}. 



The orthogonal projection on ker(6 • Vy) (with respect to the scalar product of L 2 (R m )), 
denoted by (•), reduces to average along the characteristic flow Y cf. [3] Propositions 2.2, 
2.3. 

Proposition 3.1 For any function u G L 2 (M m ) the family (u) T := ^ u(Y(s; -))ds,T > 
converges strongly in L 2 (M m ), when T — > +oo, towards the orthogonal projection of u on 
ker (6 • V y ) 

lim (u) T = (u) , (it) G ker(6 • Vy) and / (it — (u))<p dy = 0, V 92 G ker(o • Vy). 
Since • Vy is antisymmetric, one gets easily 

Range (b ■ %) = (ker(6 • Vy)) 1 = ker(Proj ker(fe . Vy) ) = ker (•) . (27) 



Remark 3.1 If u G L 2 (R m ) satisfies J Rm u(y)b • VyV> dy = 0, V ip G C£(R m ) and J Rm it<£ dy = 
0, V <p G ker(6 • Vy), f/ten it = 0. Indeed, as u G L 2 (R m ) C Lj 0C (R m ), f/te /irsi condition says 
that b ■ V y u = in D'(M m ) and thus u G ker(o • Vy). Using now the second condition with 
ip = u one gets J Rm ^ 2 dy = and i/ms u = 0. 

In the particular case when Range (b ■ Vy) is closed, which is equivalent to the Poincare 
inequality (cf. [7j pp. 29) 

3C P >0 : (( (u-{u)) 2 dy) ' < C P ( [ (b-V y u) 2 dy) ' , it G dom(6 • Vy) (28) 
(|27|) implies the solvability condition 

3 it G dom(6 • Vy) such that b ■ V y u = v iff (u) = 0. 
If || • || stands for the L 2 (R m ) norm we have 

Proposition 3.2 Under the hypothesis (|28p . b ■ Vy restricted to ker (•} is one to one map 
onto ker (•). Its inverse, denoted (b • Vy) -1 , belongs to £(ker (•} ,ker (•)) and 

||(6 • V v )~ ||,C(ker(->,ker(-» < Cp- 

Another operator which will play a crucial role is T = —div y (b <8> oVy) whose domain is 

dom(T) = {u G dom(6 • Vy) : b ■ V y u G dom(6 • V y )}. 

The operator T is self-adjoint and under the previous hypotheses, has the same kernel and 
range as b ■ Vy. 
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Proposition 3.3 Under the hypotheses (fl3j) . (fT3|) . ([28]) i/ie operator T satisfies 

ker T = ker(6 • V y ), Range T = Range (b ■ V y ) = ker (•} 
and \\u — (u) \\ < Cp\\Tu\\,u G dom(T). 

Proof. Obviously ker (6 ■ V y ) C kerT. Conversely, for any u G kerT we have J" Hm (b ■ 
V y u) 2 dy = L m uTu dy = and therefore u G ker (6 • V y ). 

Clearly Range T C Range (6 • V y ) = ker (•}. Consider now w G ker (•} = Range (6 • Vy). 
By Proposition 13.21 there is v G ker (•} ndom(6- V y ) such that b-V y v = w. Applying one more 
time Proposition [321 there is u G ker (•) n dom(6 • V y ) such that b ■ V y u = v. We deduce that 
u G domT, w = T(—u). Finally, for any u G domT we apply twice the Poincare inequality, 
taking into account that (b ■ V y u) = 

\\u - (u) || < C P \\b ■ V y u\\ < Cp\\T u\\. 

□ 

Remark 3.2 The average along the flow ofb can be defined in any Lebesgue space L g (W m ), 
q G [1, +oo]. We refer to J3J/ for a complete presentation of these results. 

3.1 Average and first order differential operators 

We are looking for first order derivations commuting with the average operator. Recall that 
the commutator • 'V y ,r] ■ V y ] between two first order differential operators is still a first 
order differential operator, whose vector field, denoted by [£, rj] , is given by the Poisson bracket 
between £ and rj 

[f ' V y ,r? • V y ] := £ ■ V y (rj ■V y )-r 1 - V a (£ • V y ) = • V y 

where [£, 77] = - V y )r] — {rj ■ V y )^. The two vector fields ^ and r\ are said in involution iff 
their Poisson bracket vanishes. 

Assume that c{y) is a smooth vector field, satisfying c(Y(s; y)) = d y Y(s; y)c(y), s G K, y G 
]R m , where Y is the flow of b (not necessarily divergence free here). Taking the derivative 
with respect to s at s = yields (6 • V y )c = d y b c(y), saying that [b, c] = 0. Actually the 
converse implication holds true and we obtain the following characterization for vector fields 
in involution, which is valid in distributions as well (see Appendix [A] for proof details). 
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Proposition 3.4 Consider b G W i '° (W n ) (not necessarily divergence free), with linear 
growth and c G L[ oc (R m ). Then (b ■ V y )c -d y bc = inV '(R m ) iff 

c(Y(s;y)) = d y Y(s;y)c(y), sel, y G M m . (29) 

We establish also weak formulations characterizing the involution between two fields, in 
distribution sense (see Appendix [A] for the proof). The notation w s stands for w o Y(s; ■). 

Proposition 3.5 Consider b G W 1 *£°(M m ), with linear growth and zero divergence and c G 

L? D (R m ). Then the following statements are equivalent 

1. 

[b,c] = mP'(R m ) 

2. 

[ {c-V y u)v- s dy= [ (c-V y u s )vdy, Vm,d £ C*(IR m ) (30) 

3. 

! c-Vyub-VyV dy+ [ c ■ V y (b ■ V y u)v dy = 0, VuGC c 2 (R m ), v G C^(IR m ). (31) 

Remark 3.3 If [b,c] = in D'(R m ), applying (|30p urai/i v = 1 on the support of u s (and 
therefore v- s = 1 on the support of u) yields 

/ c-V y udy= / c-VyU s dy, ueC^(R m ) 

saying that div y c is constant along the flow of b (in T>'(R m )). 

We claim that for vector fields c in involution with b, the derivation c • V y commutes with 
the average operator. 

Proposition 3.6 Consider a vector field c G L* oc (R m ) with bounded divergence, in involution 
with b, that is [b,c] = in T>'(R m ). Then the operators u — >• c • V y u, u — > divy(uc) commute 
with the average operator i.e., for any u G dom(c • V y ) = dom(div J/ (- c)) we have (u) G 
dom(c • Vy) = dom(div y (- c)) and 

(c • Vyti) = c • V y (u) , (divy(uc)) = d\v y ((u) c). 
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Proof. Consider u G dom(c • V y ), s£R and </? G C^M" 1 ). We have 

/ u s c • V y (/? dy = u{c ■ V y <p)- S dy (32) 

jR m Jl m 



u(c • V y )c£_ s dy 
div y (uc)(p- s dy 

I 

(div y (uc)) s ip(y) dy 



saying that u s G dom(c • V y ) = dom(divy(- c)) and div y (u s c) = (div y (uc)) s . We deduce 
c • Vy-Us = (c • V y u) s cf. Remark 13.31 Integrating (j32|) with respect to s between and T > 
one gets 



/ — / u s ds c • V y y? dy = — / I u s c - V y ip dyds 
jR ml Jo 1 Jo Jm. m 



1 

T 



o 



T 

(divj / (uc)) s (^(y) dyds 



1 

By Proposition 13.11 we know that ^ u s ds — > (u) and y J T (div y (nc)) s ds — > (div^uc)) 
strongly in L 2 (M m ), when T — > +oo, and thus we obtain 



o 



T 

(div y («c)) s ds (p(y) dy. 



(u) c ■ Vytp dy = - / (dw y (uc)} <p(y) dy 

■>■ JR™ 

saying that (u) G dom(c • V y ) and div y ((-u) c) = (divj / (uc)}, c • V y (u) = (c • V y u). 
3.2 Average and second order differential operators 

We investigate the second order differential operators — d\v y {A{y) 1 \J y ) commuting with the 
average operator along the flow of b, where A(y) is a smooth field of symmetric matrix. Such 
second order operators leave invariant ker(6- V y ). Indeed, for any u G dom(— dxv y {A(y)V y )) n 
ker(6 ■ V y ) we have 

-div y (A(y)Vyu) = -div y {A(y) (u)) = (-div y (A(y)X7 y u)) G ker(6- V y ). 

For this reason it is worth considering the operators — div y (A(y)\/ y ) commuting with b ■ Vy. 
A straightforward computation shows that 

Proposition 3.7 Consider a divergence free vector field b G W 2 ' oc (W n ) and a matrix field 
A G W 2,0O (W n ). The commutator between b ■ V y and —div y (A(y)V y ) is still a second order 
differential operator 

[b ■ Vy, -diVy(AX7y)] = ~ dW y ([b , A]V y) 
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whose matrix field, denoted by [b,A], is given by 

[b, A] = (b- V y )A - d y bA(y) - A(y) %b, y G R m . 

Remark 3.4 We have the formula 1 [b, A] = [b, t A] . In particular if A(y) is a field of symmet- 
ric (resp. anti- symmetric) matrix, the field [b,A] has also symmetric (resp. anti- symmetric) 
matrix. 

As for vector fields in involution, we have the following characterization (see Appendix lAl for 
proof details). 

Proposition 3.8 Consider b G W 1 *' ( ?°(M m ) (not necessarily divergence free) with linear growth 
and A(y) G L\ oc (W m ). Then [b,A] = in V'(R m ) iff 

A(Y(s;y)) = d y Y(s;y)A(y) t dyY(s;y), s G R, y G M m . (33) 

For fields of symmetric matrix we have the weak characterization (see Appendix [X] for the 
proof). 

Proposition 3.9 Consider b G W / 1 ^ x> (M m ) with linear growth, zero divergence and A G 

L 1 1 oc (M m ) a field of symmetric matrix. Then the following statements are equivalent 

1. 

[b,A] =0 in V'(R m ). 

2. 

/ A(y)V y u s ■ V y v s dy = A(y)V y u ■ V y v dy 
jr™ Jm. m 

for any s G M, u,v G C^M" 1 ). 
3. 

/ A{y)Vy(b-Vyu)-V y vdy+ I A(y)V y u-V y (b-V y v)dy = 
for any u,v G C^(R m ). 

We consider the (formal) adjoint of the linear operator A — > [b, A], with respect to the 
scalar product (U, V) = J„ m U(y) : V(y) dy, given by 

Q^-(b- V y )Q - 'dybQiy) - Q(y)d y b 

when diVyb = 0. The following characterization comes easily and the proof is left to the 
reader. 
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Proposition 3.10 Consider b G W l '°°(R m ), with linear growth and Q G L^iW 71 ). Then 
-(b-Vy)Q- t d y bQ(y)-Q(y)d y b = mV'{R m ) iff 



Remark 3.5 If Q(y) satisfies ([HI"]) and is invertible for any y G R m with Q 1 G Lj^IR" 1 ), 
f/ien = d v Y(s;y)Q- 1 (y) t d y Y(s;y), s G R,y G R m and therefore [6, Q -1 ] = in 

D'(R m ). If P(y) satisfies (j33|) and is invertible for any y G R m ; £/zen 

p-^y^y)) = t ^ 1 ( S ;y)p- 1 (y)aj / y- 1 ( S ;y), set, y G R m 

and therefore -(b ■ V y )P - t d y bP{y) - P(y)d y b = inV'(W n ). 

As for vector fields in involution, the matrix fields in involution with b generate second order 
differential operators commuting with the average operator. 

Proposition 3.11 Consider a matrix field A G L\ oc (R m ) such that div y ^4 G L\ oc (R m ) and 
[b,A] = in T>'(R m ). Therefore the operator u —> — &iv y (AV y u) commutes with the average 
operator i.e., for any u G dom(— dxv y (AVy)) we have (u) G dom(— div y (A'V y )) and 

- (divy(AVyu)) = -div y (AVy (u)). 

Proof. Consider u G dom(-div y (^V y )) = {w G L 2 (R m ) : -div y {AV y w) G L 2 (R m )}. For 
any s G R, <p G C c 2 (R m ) we have 



Q(Y(s;y)) 



l dyY \s;y)Q(y)d y Y- 1 (s;y), set, y G R' 



(34) 




(35) 



By the implication 1. 



2. of Proposition 13.91 (which does not require the symmetry of 



A{y)) we know that 




for any ij> G C c 2 (R m ). We deduce that 




and thus (div 3/ ( 



^AVytp^s = div y (fAVy(p- s 



). Combining with ([35]) yields 




(36) 



- J Rm div y (AVyu)ip- s dy 
f Rm (div y (AV y u)) s cp(y) dy 
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saying that u s £ dom(— div y (A S/ y )) and 

-div J/ (ylVj / 'u s ) = (-divy(AV y u)) s . 
Integrating (f36j) with respect to s between and T we obtain 

/ — / u s ds div^f^V^) dy = I ^ [ (dw y (AV y u)) s ds <p(y) dy. 

JWl™ 1 JO JR ml JO 

Letting T — > +00 yields 



and therefore (u) G dom(div 3/ (^4V 3/ )), div y (AS/ y {u}) = {div y (AV y u)) . |— j 
3.3 The averaged diffusion matrix field 

We are looking for the limit, when e — > 0, of {TJ, ([2]). We expect that the limit u = lim^o u £ 
satisfies (0), (0). By (0) we deduce that at any time t € R + , •) € ker(6- V y ). Observe also 
that div y (6(8>6Vj / n 1 ) = b-V y {b-V y u l ) 6 Range {b-V y ) C ker (•) and therefore the closure for 
u comes by applying the average operator to ([6|) and by noticing that {dtu) = dt (u) = dtu 



At least when [b, D] = 0, we know by Proposition 13.111 that 

{divy(DV y u)) = dxv y {DV y {u)) = d\v y {DV y u) 

and (f37|) reduces to the diffusion equation associated to the matrix field D(y). Nevertheless, 
even if [b, D] ^ 0, ()37p behaves like a diffusion equation. More exactly the L 2 (]R m ) norm of 
the solution decreases with a rate proportional to the L 2 (M m ) norm of its gradient under the 
hypothesis ([3]) 




d t u - (dw y {DV y u)) =0, ye K ; 



(37) 
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We expect that, under appropriate hypotheses, (f37|) coincides with a diffusion equation, 
corresponding to some averaged matrix field V, that is 

3V(y) : [b,T>]=0 and (-div„(DV w u)) = -div y (X>V„«), Vm£ ker(6 • V y ). (38) 

It is easily seen that in this case the limit model (j37|) reduces to 

d t u - div„(DV v «) = 0, i G R+, t/el m . 

In this section we identify sufficient conditions which guarantee the existence of the matrix 
field T>. We will see that it appears as the long time limit of the solution of another parabolic 
type problem, whose initial data is D, and thus as the orthogonal projection of the field D(y) 
(with respect to some scalar product to be defined) on a subset of {^4 G L 1 1 oc (lR m ) : [b,A] = 
in ©'(R" 1 )}. We assume that (fTT|) holds true. We introduce the set 

H Q = {A = A{y) : / Q(y)A{y) : A(y)Q(y) dy < +00} 

where Q = P^ 1 and the bilinear application 

(., -)q : H Q x H Q R, (A,B) Q = I Q{y)A(y) : B(y)Q{y) dy 

JR™ 

which is symmetric and positive definite. Indeed, for any A G Hq we have 

(A,A) Q = [ Q 1/2 AQ 1/2 : Q 1/2 AQ 1/2 dy > 

with equality iff Q 1 ! 2 AQ 1 ! 2 = and thus iff A = 0. The set Hq endowed with the scalar 

1/2 

product (•, -)q becomes a Hilbert space, whose norm is denoted by \A\q = (A, A)a , A G Hq. 
Observe that Hq C {A(y) : A € Lj L oc (]R m )}. Indeed, if for any matrix M the notation \M\ 
stands for the norm subordonated to the euclidian norm of R m 



\ M \= sup ^TTp- < (M : M) 1 ^ 2 
^\{o} Isl 



we have for a. a. y G 



1/1/ M A (y)€-V ,o n \ 

\A(y)\ = sup (39) 

£,^0 If I M 

Ql/2 A Ql/2pl/2^ . pl/2 lpl/2^1 [pl/2 | 
sup " ' 



< |Q 1/2 AQ 1/2 | |pV2|2 

< (g 1 / 2 ^ 1 / 2 : qWaq 1 /*) 1 ' 2 \P\. 

We deduce that for any i? > 

/ \A(y)\ dy < [ (Q 1/2 AQ 1/2 : Q^AQ^f/ 2 \P\ dy < (A, A) 1 ^ ( [ \P(y)\ 2 dy 
Jb r Jb r \Jb r j 



1/2 
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Remark 3.6 We know by Remark \3.5\ that Q s = t d y Y~ 1 (s;y)Q(y)d y Y~ 1 (s;y) which writes 
t O(s;y)0(s;y) = I where 0(s;y) = Q^J 2 d y Y(s;y)Q~ l l 2 . Therefore the matrix 0(s;y) are 
orthogonal and we have 

Q 1 s / 2 d y Y(s;y)Q- 1 / 2 = 0(s;y) = ^(ajy) = Q~ 1/2 'dyY^Q 1 ' 2 (40) 

Q- 1 ' 2 t d y Y(s;y)Q 1 J 2 = = 0-\s;y) = Q l/2 d y Y~ 1 Q~ 1 / 2 . (41) 

Proposition 3.12 The family of applications A — > G(s)A := d y Y~ 1 (s; -)A S t d y Y~ 1 (s; ■) is 
a C°- group of unitary operators on Hq. 

Proof. For any A G Hq observe, thanks to (|4ip . that 
\d y Y-\s;-)A s t d y Y- 1 (s;-)\i = [ Q l/2 d y Y~ 1 A s t d y Y~ 1 Q l/2 : Q^dyY' 1 A S %Y~ X Q X I 2 dy 

= 1*0(8; y)Q 1 J 2 A s Q 1 J 2 0(s; y) y)Q 1 J 2 A s Ql/ 2 0(s; y) dy 

= [ Q l J 2 A s Q\/ 2 : Q\/ 2 A S Q\/ 2 dy 
= [ Q l / 2 AQ 1 / 2 : dy 
= \A\ 2 Q - 

Clearly G(0)A = A, A £ Hq and for any s, t G R we have 

G(a)G(t)A = ^y- 1 ^; O^^A),'^^ 1 ^; •) 

= fyY" 1 ^; •)( ^ (a; ■))(A t ) s t (d y Y)- 1 (t; Y(s; -^dyY-^s; •) 
= d y Y-\t + a; O^t+Z^y- 1 ^ + s; •) = G(t + s)A, A £ Hq. 

It remains to check the continuity of the group, i.e., lim s _^o G(s)A = A strongly in Hq for 
any A £ Hq. For any s £ R we have 

|G( S )^ - A\ 2 Q = \G(s)A\ 2 Q + |^ - 2(G(s)A, A) Q = 2\A\ 2 Q - 2(G(s)A, A) Q 

and thus it is enough to prove that lim^o G(s)A = A weakly in Hq. As |G(s)| = 1 for any 
s G R, we are done if we prove that lim s ^ (G(s)A, U)q = (A, U)q for any U G C°(R m ) C Hq. 
But it is easily seen that lim s _>.o G(—s)U = U strongly in Hq, for U G C°(R m ) and thus 

lun(G(s)A, U) Q = limM, G(-s)U) Q = (A, U) Q , U G C c °(R m ). 

□ 
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We denote by L the infinitesimal generator of the group G 

G(s)A — A 

L : dom(L) C Hq ->■ Hq, domL = {A £ Hq : 3 lim — — in Hq} 

s— >0 S 

and L(A) = lim s _> G(s) s A ~ A for any A G dom(L). Notice that C c 1 (M m ) C dom(L) and L(A) = 
b ■ V y A - dybA - A t d y b, A G C c 1 (M m ) (use the hypothesis Q G Lf oc (M. m ) and the dominated 
convergence theorem). Observe also that the group G commutes with transposition i.e. 
G(s) l A = t G(s)A, s G Hq and for any A G dom(L) we have l A G dom(L), L^A) = 

t L{A). The main properties of the operator L are summarized below (when b is divergence 
free). 

Proposition 3.13 

1. The domain of L is dense in Hq and L is closed. 

2. The matrix field A G Hq belongs to dom(L) iff there is a constant C > such that 

\G{s)A-A\ Q < C\s\, s£R. (42) 

3. The operator L is skew-adjoint. 
4- For any A G dom(L) we have 

-div y (L(A)Vy) = b ■ Vy(-div y (AVy)) + dh y (AV y (b ■ V„)) in V '(R m ) 

that is 

/ L(A)V y u ■ V y v dy = - AV y u ■ V y (b ■ V y v) dy - AV y (b ■ V y u) ■ V y v dy 



for any u,v G C c 2 (R m ). 

Proof. 1. The operator L is the infinitesimal generator of a C°-group, and therefore dom(L) 
is dense and L is closed. 

2. Assume that A G dom(L). We know that 4;G(s)A = L(G(s)A) = G{s)L(A) and thus 



\G(s)A-A\ 



G{t)L{A) dr 



o 



< 







\G(t)L(A)\q dr 



\s\ \L(A)\ Q , sG 



Conversely, assume that (j4"2j) holds true. Therefore we can extract a sequence (sk)k converging 
to such that 

lim *~*^ Sfc ^ = V weakly in Hq. 

k— »+oo Sk 

For any U G dom(L) we obtain 

(G(s k )A -A \ = / G(-s k )U - U 



\ s k J Q V Sk 
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and thus, letting k — y +00 yields 



(V,U) Q = -(A,L(U))c 



(43) 



But since U G dom(L), all the trajectory {G(t)U : r G M} is contained in dom(L) and 
G{-s k )U = U + J Q Sh L(G{r)U)dT. We deduce 



(G(s k )A - A, U) Q = L(G(r)U) dr 



"Sfc 



(A,L(G(t)U)) q dr 



Taking into account that 



/- Sfc G(r)Udr 



G(s k )A - A 



U 



Sk 



S \v,G(t)U) q dr 



V,J " G(t)U dr 



^ l s fcl \U\q we obtain 



< \V\q\U\q, U G dom(L) 



and thus, by the density of dom(L) in Hq one gets 

G(s k )A-A 



Sk 



<\v\ Q , ken. 



Since V is the weak limit in Hq of ( G ^ Sk }^ — ) , we deduce that lim^ +00 ^^ k >^ - - ( 



G(s k )A-A 



strongly in Hq. As the limit V is uniquely determined by ()43|) . all the family [ QislAzA 
converges strongly , when s — )■ 0, towards V in Hq and thus A £ dom(L). 
3. For any U, V € dom(L) we can write 

(G(s)C/ - C/; F)q + (U,V- G(-s)V)q = 0, s € R. 

Taking into account that 

lhn mu-u = ihn v- a{ - s)v = 

s^O S s^O S 

we obtain (L(U),V) Q + (U,L(V))q = saying that V G dom(L*) and L*(V) = -L(V). 
Therefore L C (— L*). It remains to establish the converse inclusion. Let V G dom(L*), i.e., 
3C > such that 

|(L(C/),y) Q | < C\U\ Q , Uedom(L). 
For any s G R, ?7 G dom(L) we have 

( G ( s )y _ V , U) Q = (V, G(-s)U - U) Q = (V, I ' LG(t)U dr) Q = / '(V, LG{t)U)q dr 

Jo Jo 
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implying 

\(G(s)V -V,U) Q \ < C\s\ \U\ Q , s GE. 

Therefore \G(s)V — V\q < C|s|,s G E and by the previous statement V G dom(L). Finally 
dom(L) = dom(L*) and L*(V) = —L(V), V G dom(L) = dom(L*). 
4. As L is skew-adjoint, we obtain 

- / L(A)V y u ■ V y v dy = —(L(A), Q^VyV ® V^nQ- 1 ) Q = (A, L(Q^ 1 V y v ® V^uQ" 1 ) ) Q . 

Recall that P = Q" 1 satisfies L(P) = 0, that is, G(s)P = P, s G E and thus 

r/ ^_i„ „ ^ n ,. G(s)PV v v ® V v uP - PV„v (8) V v uP 
L(Q _1 Vj,v<8>Vj,itQ _1 ) = lim — ^ * * * y - — 



lim 



s^O S 

d y Y- 1 (s- ■)P s (Vyv) s <8 (V^^.P^y- 1 ^; •) - PV^ ® V„iiP 



S 

_ P l dyY{s; ■){V y v) s <g> (V y u) s d y Y(s; -)P - PV y v g V^P 

s-s>0 S 

- lim PV ^ S ® V 2/ UsP ~ PV 2/' U ® V^uP 

s-5>0 S 

= PVy(6 ■ VyV) <g> VyUP + PV yV <g> V„(& ■ V yU) P. 

Finally one gets 

- / L(A)V y u ■ V y v dy = (A, PV y (b ■ V y v) <g> V y uP) + PV y v <g> V y (b ■ V y u)P) Q 

= AVyU- V y (b- V y v) dy+ / AV y (b ■ V y u) ■ V y v dy. 

□ 

We claim that dom(L) is left invariant by some special (weighted with respect to the ma- 
trix field Q) positive/negative part functions. The notations A ± stand for the usual posi- 
tive/negative parts of a symmetric matrix A 

A ± = SA ±t S, A = SA t S 

where A, A 1 * 1 are the diagonal matrix containing the eigenvalues of A and the positive/negative 
parts of these eigenvalues respectively, and S is the orthogonal matrix whose columns contain 
a orthonormal basis of eigenvectors for A. Notice that 

A + : A' = 0, A + - A' = A, A + : A + + A~ : A~ = A : A. 

We introduce also the positive/negative part functions which associate to any field of sym- 
metric matrix A(y) the fields of symmetric matrix A ( ^ ± (y) given by 

Ql/2 A Q± Ql/2 = (g l/2 Ag l/2 )± _ 

Observe that A Q+ - A Q ~ = A. 
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Proposition 3.14 

1. The applications A — > leave invariant the subset {A € dom(L) : 1 A = A}. 

2. For any A £ dom(L),*j4 = A we have 

{a Q+ jA Q~ )q = 0) ( L ( A Q+) 5 L{AQ-)) Q < 0. 
Proof. 1. Consider A G dom(L),*j4 = A. It is easily seean that t A ( ^ ± = A^ and 
\A®% + \A<3-\1 = I (Q 1/2 AQ 1/2 ) + : (QV*AQV 2 ) + dy 

jR m 

+ / (Q^AQ 1 / 2 )- : (Q^AQ 1 / 2 )- dy 

JR m 

= / Q 1/2 AQ 1/2 : Q 1/2 AQ 1/2 dy = \A\ 2 Q < +oo 

JR m 

and therefore A^ € i^Q. The positive/negative parts A^ are orthogonal in Hq 
( A Q+,aQ-) q = [ (Q^AQ^Y ■■ (Q 1/2 AQ 1/2 r dy = 0. 

JR m 

We claim that A^ satisfies (|42|) . Indeed, thanks to (|41|) we can write, using the notation 
X- 2 = X : X 

| G ( s)j4 Q± _ A Q±\2 Q = f {Q^idyY-'iA^/dyY- 1 - A Q± )Q 1/2 } :2 dy (44) 

jR m 

= f { t O{s-y)Q 1 J 2 (A Q± ) s Q 1 J 2 0{s;y) - Q 1/2 A Q± Q 1/2 } :2 dy 

JR m 

= f {^{s-y^AsQll^Ois^-^AQy^rdy. 

JR m 

Similarly we obtain 

\G(s)A - A\ 2 Q = [ {*0( S ; y)Ql/ 2 A s Q l J 2 0(s; y) - Q 1/2 AQ 1/2 }'- 2 dy. (45) 

JR m 

We are done if we prove that for any symmetric matrix U, V and any orthogonal matrix R 
we have the inequality 

( t RU ± R - V ± ) : ( t RU ± R -V ± )<( l RUR - V ) : ( t BUR -V). (46) 

For the sake of the presentation, we consider the case of positive parts U + ,V + . The other 
one comes in a similar way. The above inequality reduces to 

2 t RUR : V - 2 t RU + R : V + < t RU^R : t RU^R + V~ : V~ 

or equivalently, replacing U by U + — U~ and 1/ by V + — V~, to 

-2 ^CZ+i? : V" - 2 : y+ + 2 *l«7~ii : V~ < t RU~R : t RU~R + V~ : V"*. 
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It is easily seen that the previous inequality holds true, since t RU + R : V > 0, RU R 
V + > and 

2 tRU-R : V~ < 2( t RU-R : l RU- R) 1/2 (V- : V') 1 / 2 < l RU-R : t RU-R + V~ : V~ . 
Combining (]3I|) . (05]) and (06} with 

t/ = Qy 2 ^ s Qy 2 , V = Q 1 ' 2 AQ 1 / 2 , R = 

yields 

sup J-AJ ^ < sup 1 y ' ^ < |L(A)| 

saying that A Q± G dom(L). 

2. For any A G dom(L), '4 = A we can write 

(A«+,A«-)q = / QWAQ+QV 2 : Q^A^-Q 1 / 2 dy 
Jm. m 

= [ (Q l / 2 AQ l / 2 ) + : (Q 1/2 AQ 1/2 y dy = 0. 

Since A^ G dom(L) we have 

s^O s 

and therefore, thanks to (01]), we obtain 



(L(AQ+),L(AQ-)) q = lim 



V / Q 

/■ Q 1/2 ( -G(-f)^ Q + )Q 1 / 2 Q x / 2 ( )g x / 2 
lim / : dy 

s^O J Rm S S 

t O(|;y)(Qf A f Qf )+0(f;y)-*O(-f;y)(Qi / jA_ f Qi / |)+O(-f;y) 

lim 

s->0 

*0(|; y)(Ql /2 A f Qi /2 )-0(f ;y) - *0(-§; y)(Q^-f Q^)"0(-f ; y) 

2 j 2 2 j 2 _ 

s 



dy 



lim / 2 ±—2 5 ? i 2 - dy 



*o(-f;y)(Q^-#Q v I) + o(-|;y):'o(f;y)(Qi /2 ^Qi /2 )-o(|;y) 



- lim / = a a ^ ^ — dy 

< 



since 



J 0(± s /2; .)(Q 1/2 ^Q 1/2 )i, /2 0(±V2; •) > 0, '0(^/2; •)(Q 1/ W /2 )^ /2 0(^/2; •) > 0. 



□ 
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We intend to solve the problem ()18p . (|19p by using variational methods. We introduce the 
space Vq = dom(L) C Hq endowed with the scalar product 

((A, B)) Q = (A, B) Q + (L(A), L(B))q, A, Be Vq. 

Clearly (Vq, ((•, -))q) i s a Hilbert space (use the fact that L is closed) and the inclusion 
Vq C Hq is continuous, with dense image. The notation || • \\q stands for the norm associated 
to the scalar product ((•, -))q 

\\Af Q = ((A,A)) Q = (A,A)q + (L(A),L(A))q = \A\ Q + \L(A)\ 2 Q , A € Vq. 

We introduce the bilinear form a : Vq x Vq — )• R 

(t(AS) = (L(A),L(S))q, A^GVq. 

Notice that a is coercive on Vq with respect to Hq 

a(A,A) + \A\ 2 Q = \\Af Q , A € Vq. 

By Theorems 1,2 pp. 620 [8] we deduce that for any D G Hq there is a unique variational 
solution for (USD, (USD that is A £ C b (R+; Hq) n £ 2 (R+; Vq), d t A G L 2 (R+; Vq) 

A(0)=A ^(A(t),U) Q +a(A(t),U)=0, in D'(R m ), V H £ Vq, 

The long time limit of the solution of f)18[) . (|19p provides the averaged matrix field in fj38[) . 
Proof, (of Theorem 12. ip The identity 

~\A(t)\% + \L(A(t))\h = 0, ^R + 

gives the estimates 

r+oo i 

|^(t)lo < PIq, *eR + , y o |L(A(t))| Q dt<-|D| Q . 

Consider (£&)& such that t k — >• +oo as — )• +oo and (A(tk))k converges weakly towards some 
matrix field X in Hq. For any U £ kerL we have 

A(^),C/)q=0, tel + 

and therefore 

(Proj kcrL AC/)Q = (AtOg = (M0),U) Q = (A(t k ),U) Q = (X,U) Q , U £ kerL. (47) 
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Since L(A) G L 2 (R + ; Hq) we deduce that lirm 3 _ > + 00 L(A(tk)) = strongly in Hq. For any 
V G Vq we have 

(X,L(F))q= lim (A(t fc ),L(^)) = - Urn (L(A(i fc )), F) Q = 0. 

k—¥+O0 K— > + 0O 

We deduce that X € dom(L*) = dom(L) and L(X) = 0, which combined with ()47|) says that 
X = Proj kerL D, or X = (D)q. By the uniqueness of the limit we obtain limj^+oo A(t) = 
Proj ker £-D weakly in Hq. Assume now that t D = D. As L commutes with transposition, we 
have dt A — L(L(*A)) = 0, 'A(O) = D. By the uniqueness we obtain * A = A and thus 

t (D) Q = '(w- lim A(i))=w- lim *A(i) = w - lim A(t) = (-D)q . 

Suppose that D > and let us check that (-D)q > 0. By Proposition 13.141 we know that 
A ( 3 ± (t) £ Vq, t G R + and 

(A«+(t) J AO-(t)) o = 0, (L(A«+(t)),L(AQ-(t))) Q <0, i € R+. 

It is sufficient to consider the case of smooth solutions. Multiplying (fTSj) by — A®~(t) one 
gets 

IA|^Q- (t) |2 + \L{A®-{t)\ 2 Q = (d t AQ+,A<3-(t)) Q + (L(A«+(t)),L(A^))) Q (48) 

<{d t A^,A^{t)) Q . 

But for any < ft < i we have 

_ A«+(t - ft), A«-(t)) g = -(A«+(t - h),A®-(t)) Q < 

and therefore (d t A Q+ (t),A Q ~(t))Q < 0. Observe that Q 1 / 2 A Q ~(0)Q 1 / 2 = {Q 1 / 2 DQ 1 / 2 )- = 
0, since Q 1 ' 2 DQ 1 / 2 is symmetric and positive. Thus A^~(0) = 0, and from (|48|) we obtain 

\\A Q -{t)\i<\\A®-m Q = v 

implying that Q 1 / 2 A(t)Q 1 / 2 > and A(t) > 0, t G R+. Take now any U & H Q , *!/ = U, 
U > 0. By weak convergence we have 

({D)q,U) q = lim (A(i),[/) g = hm / Q 1/2 A(f)Q 1/2 : Q 1/2 UQ 1/2 dy > 

and thus (D)q > 0. By construction (D)q = Proj keri -D G kerL. It remains to justify the 
second statement in (|2(jp . and (|21|) . Take a bounded function ip G L°°(R m ) which remains 
constant along the flow of b, that is (p s = if, s G R, and a smooth function u G C 1 (W m ) such 
that u s =«,s£l and 

/ (V y u ■ Q^Vyu) 2 dy < +oo. 
23 



We introduce the matrix field U given by 

U(y) = <p{y)Q-Hv) V„« ® V„« y € R m . 

By one hand notice that U G Hq 

\U\ 2 Q = [ Q 1 ' 2 UQ 1 l 2 :Q 1 l 2 UQ 1 l 2 dy = [ p^Q'^V y u\Uy 

< \\<P\\l°° I {VyU-Q^Vyuf dy. 
Jm. m 

By the other hand, we claim that U G kerL. Indeed, for any selwe have 

V y u = V y u s = t d y Y(s; y)(V y u) s 

and thus 

Q S U S Q S = <p s (VyU) s <g> (VyU) s 

= CdyY^VyU) ® COyY^VyU) 

= If l dyY~ l VyU ® VyU dyY' 1 

= t d y Y' 1 QUQdyY- 1 . 

Taking into account that Q s = t d y Y~ 1 Qd y Y~ l we obtain 

t dyY- 1 Qd„Y- 1 Us t d„Y- 1 QdyY- 1 = tdyY^QUQdyY- 1 

saying that U s (y) = d y Y(s;y)U(y) t d y Y(s;y). As (D) Q = Proj kcrL Z) one gets 

= (D - (D) Q , U) Q = [ (D- (D) Q ) : QUQ dy 

= [ <p(y)(D-(D) Q ):Vy U ®Vyudy 

= / ip(y){V y u ■ DVyU - V y u ■ (D) Q V y u} dy. 

In particular, taking ip = 1 we deduce that V y u ■ (D)q V y u G L 1 (lR m ) and 

/ V y u ■ (D) Q V y u dy = / V y u ■ DVyU dy = (D^' 1 V y u® V y u Q~ 1 ) Q < +oo 

since D G Hq, Q~ 1 V y u <S> V y uQ~ l G Hq. Since (D)q G kerL, the function V y u ■ (D)q V y u 
remains constant along the flow of b 

(VyU) S ■ ({D) Q ) s (VyU) S = (VyU) s ■ dyY fry) (£>) Q ' 8yY (s ' J/) (V yU) , = Vy U -(D) Q VyU. 

Therefore the function V y u ■ (D)q V y u verifies the variational formulation 

VyU ■ (D)q VyU G L 1 (M m ), (VyU ■ (D) Q VyU) s = VyU ■ (D) Q VyU, s£R (49) 
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and 



/ V y u- DV y uip dy = V y u- (D) Q V y ULpdy, M p G L°°(R m ), p s = ip, s el (50) 

It is easily seen, thanks to the hypothesis D G L°°(R m ), that ()49|) . (|50p also make sense for 
functions u G -ff^R*™) such that u s = «, s £ I. We obtain 

VyU • (D) Q V y u = (V y u ■ DVyu) , u G il^R™), u s = u, seR 

where the average operator in the right hand side should be understood in the L 1 (R m ) setting 
cf. Remark IO Moreover, i£u,v€ H 1 ^™) nker(6 • V y ) then <D)q /2 V y u, (£>}q /2 V tf u belong 
to L 2 (R m ) implying that V y u-(D) Q V y v G L 1 (W m ). As before we check that V y u-(D) Q V y v 
remains constant along the flow of b and for any ip G L°°(R m ), ip s = ip, s G R we can write 

2 / VyU ■ DVyV tp dy = / Vj,(ii + v) • DV y (u + v) ip dy 

— / V y u ■ DV y u ip dy — V y v ■ DVyV (p dy 
= / V y (u + v) ■ (D) Q V y (u + v) tp dy 

- / VyU ■ (D) Q VyU ipdy - V y v ■ (D) Q V y v ip dy 
= 2 VyU ■ (D) Q V y v ip dy. 

Finally one gets 

VyU ■ {D)q V y v = (VyU-DVyv), u , v G H 1 (K" ) n ker (6 • V y ) . 

Consider now u G F 1 (R m )nker(6-V y ) and ip G C 2 (R m ). In order to prove that (v y u ■ {D) Q V y (b ■ V y ip 
0, where the average is understood in the L 1 (R m ) setting, we need to check that 

f ip(y)VyU-(D) Q Vy(b-V y iP)dy = 

JR m 

for any ip G L°°(R" 1 ), tp s = ip, s G R. Clearly B(y) := ip(y) (D)q (y) G ker L and therefore it 
is enough to prove that 

/ V y u ■ BV y (b ■ Vyip) dy = 
for any B G kerL, which comes by the third statement of Proposition 13.91 | — | 



Remark 3.7 Assume that there is uq satisfying uo(Y(s;y)) = u${y) + s, s G R, y G R m . 
Notice that uq could be multi-valued function (think to angular coordinates) but its gradient 
satisfies for a. a. y G R m and s G R 

V y u = t d y Y(s;y)(VyU )s 
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exactly as any function u which remains constant along the flow of b. For this reason, the 



en 



last equality in (|20p holds true for any u, v £ H (M m ) D ker(6 • V„) U {no}- In the case wh 
m — 1 independent prime integrals of b are known i.e., 3u\, ...,u m -i £ H 1 (W n ) nker(6 ■ Vy), 
the average of the matrix field D comes by imposing 

VyUi ■ (D)q VyUj = (WyUi ■ DVyUj) , l,j £ {O, ...,171 ~ 1} . 

4 First order approximation 

We assume that the fields D(y), b(y) are bounded on W 71 

D £ L°°(R m ), &£L°°(R m ). (51) 

We solve ([I]), © by using variational methods. We consider the Hilbert spaces V := 
H 1 (M. m ) c H := L 2 (IR m ) (the injection V C H being continuous, with dense image) and 
the bilinear forms a £ : V x V —¥ R given by 

a £ (u,v)= I D(y)V y u -V y v dy + - [ (b ■ V y u) (b ■ V y v) dy, u,v£V. 

JR m ' ' e Jr™ 

Notice that for any < e < 1 and v £ V we have 

a £ (v, v) + d\v\ 2 H > I D{y)V y v ■ V y v + (b ■ V y v) (b ■ V y v) dy + d (v(y)) 2 dy 

jR m jR m 

>d I \V y v\ 2 dy + d f (v(y)) 2 dy 

JR m JR™ 

= d\ V \ 2 y 

saying that a 6 is coercive on V with respect to H. By Theorems 1,2 pp. 620 |8j we deduce that 
for any uf n £ H, there is a unique variational solution for ([!]), fl2J), that is u £ £ Cb(R + ; H) n 
L 2 {R + ;V) and 

u £ (0) = u £ n , — I u £ (t,y)v(y) dy + a £ (u £ (t) , v) = 0, mV'(R m ), V v £ V. 



dt ... 

By standard arguments one gets 

Proposition 4.1 T/ie solutions (u £ ) e satisfy the estimates 

r+oo r 

\\u £ \\ Cb (R + ;H) <\ufn\H, / / \V y U £ \ 2 dydt < 

.In ./lB«n 



and 

/ e \ 1/2 
life - v XIIl2(r +; h) < ( 2(1 - g ) J \ u ^\h, £6(0,1). 

We are ready to prove the convergence of the family (u £ ) e , when e \ 0, towards the solution 
of the heat equation associated to the averaged diffusion matrix field (D)n- 
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Proof, (of Theorem 12, 2p Based on the uniform estimates in Proposition 14, 1[ there is a 
sequence (Ek)k, converging to 0, such that 

u £k -± u weakly * in L°°(R + ; H), V y u £k V y u weakly in L 2 {R + ;H). 

Using the weak formulation of (HJ with test functions rj(t)ip(y), 7] G C^(R + ),ip G C^(W m ) 



yields 



V '{t)<p{y)u £k (t, y) dydt - 7/(0) / pug d V + / / V^yU E " ' DV y <p dydt 

Jo Jm.™ 



r](t)(b-V y u £k )(b- V y <p)dydt. (52) 



1 

e k Jo 

Multiplying by and letting k +oo, it is easily seen that 

i>+oo r 

/ / 77(6 - V y u) (b ■ Vyip) dydt = 0. 

JO JR m 

Therefore u(t, •) G kerT = ker(6 • V y ), t G M+, cf. Proposition 13,31 Clearly ([52]) holds true 
for any <p G V. In particular, for any ip D ker(6 • Vy) one gets 

- / r]'u £k tpdydt-r](0) u £k ip dy +/ / r/V y n £fe • Z)V y ^ dydt = 0. (53) 

Thanks to the average properties we have 

uf*<pdy= / p dy ->■ / u in ip dy 



and thus, letting k — >■ +00 in (|53j) . leads to 



/■ r+00 /• 

rj'uip dydt — 77(0) / iti n </? dy + / r]V y u ■ D\7 y ip dydt = 0. (54) 
Since u(t, ■), ip G V D ker(6 • V y ) we have cf. Theorem 12,11 

/ Vy-u • DV^y? dy = V y u ■ (D) Q V y ip dy 

Jr™ jR m 



and (f54"|) becomes 

hOO/> p r + OO 



-I j rj'uip dydt - rj(0) / u in ip dy + / rjSJ y u ■ (D) Q V y (p dydt = 0. (55) 
Jo Jr.™ Jr™ Jo Jr™ 

But ([55]) is still valid for test functions ip = b ■ V y ijj, ip G C% (R m ) since •) G ker(6 • V y ), 
u m = w — lim^o (uf n ) G ker(6 • V y ) and (-D)q G kerL 

/ u(t,y)b- V y ijdy = 0, / -u in 6 • dy = 0, / V y u • (D) Q V y (b ■ V y ip) dy = 

jR m Jw m Jr™ 

cf. Theorem 12 ,11 Therefore, for any v G V one gets 

£ / u(t,y)u(y)di/+ / V y u-(D) Q V y vdy = mV'(R m ) 

with it(0) = Ui n . By the uniqueness of the solution of (f22]h (j23|) we deduce that all the family 
(u £ ) e converges weakly to u. 
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Remark 4.1 Notice that ([22]) propagates the constraint b ■ V y u = 0, if satisfied initially. 
Indeed, for any v G C^(M. m ) we Ziawe 

— / u{t,y)v{y)dy + I V y u ■ (D) Q V y v dy = m V >(R m ). (56) 

Since (D)q G kerL, we fcnoro i/ie second statement of Proposition HOI tfttrf 

/ V y n s • (D) Q VyV dy= V y u ■ (D)q V y v- S dy. 

Replacing v by V- s in (|56p we obtain 

£ / u s vdy+ [ VyU s -(D) Q V y vdy = inV'(R m ) 

and therefore u s solves 

d t u s - div y ((D) Q Vyu s ) = 0, (t, y) G K+ x K m 

and u s (0, y) = ni n (y(s;y)) = u m {y),y G M m . -By the uniqueness of the solution of (|22p . 
one gets u s = u and thus, at any time t G R+, 6 • V y u{t, •) = 0. 

5 Second order approximation 

For the moment we have determined the model satisfied by the dominant term in the expan- 
sion @. We focus now on second order approximation, that is, a model which takes into 
account the first order correction term eu 1 . Up to now we have used the equations ©, ([6]). 
Finding a closure for u + eu 1 will require one more equation 

d t u l - d\-v y {DVyU 1 ) - div y (b <g> bV y u 2 ) = 0, (t, y) G K+ x R m . (57) 

Let us see, at least formally, how to get a second order approximation for {u £ ) £ , when e 
becomes small. The first order approximation i.e., the closure for u, has been obtained by 
averaging ^ and by taking into account that u G ker(6 ■ V y ) 

d t u = {divy(DVyu)) = div y {(D) Q V y u). 

Thus u 1 satisfies 

div y ({D) Q V y u) - dwy(DVyu) - div y (b <g> bVyu 1 ) = (58) 
from which we expect to express u 1 , up to a function in ker(6 • V y ), in terms of u. 
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Proof, (of Theorem 12, 3p We claim that Range L 2 = Range L and thus Range L 2 is closed 
as well. Clearly Range L 2 C Range L. Consider now Z = L(Y) for some Y G dom(L). But 



Y — Proj kcri Y G kerL = (kerL*) = Range L = Range L and there is X G dom(L) such 
that Y - Proj kerL y = L(X). Finally X G dom(L 2 ) and 

Z = L(Y) = L(Y - Proj kcrL y) = L(L(X)). 



By construction we have D — {D)q € (kerL) = (kerL*) = Range L = Range L = 
Range L 2 and thus there is a unique F £ dom(L 2 ) n (ker L) 1 - such that D = {D)q — L(L(F)). 
As F e (kerL)- 1 , there is C G dom(L) such that F = L(C) implying that t F = t L{C) = 
L(*C). Therefore t F G dom(L 2 ) n (kerL)- 1 - and satisfies the same equation as F 

L(L( t F))= t L(L(F)) = (D) Q -D. 

By the uniqueness we deduce that F is a field of symmetric matrix. By Proposition 13.131 we 
know that 

-div y (L{F)V y ) = [b- Vy,-dw y {FVy)} in V'(R m ) 

i.e., 

i L(F)V y u ■ V y v dy = - / FV y u ■ X7 y (b ■ V y v) dy - FV y (b ■ V y u) ■ V y v dy 

jR m JR™ il m 

for any u, v G C 2 (K m ). Similarly, E := L(F) satisfies 

-dwy(L 2 (F)V y ) = -dw y {L(E)Vy) = [b ■ V y , -div y (EVy)] in V'(R m ) 

and thus, for any u, v G C^(lR m ) one gets 

/ {{D)n — D)VyU ■ V y v dy = [ L 2 (F)V y u ■ V y v dy 

= - / L(F)V y u ■ V y (b ■ Vyv) dy- L(F)V y (b ■ V y u) ■ V y v dy 

jR m ' ' JR.™ 

= FV y u-Vy(b-Vy{b-Vyv)) dy+ FVy{b-Vyu)-Vy{b-V y v)dy 

JR m JR m 

+ FVy{b-V y u)-Vy{b-V y v)dy+ FV y (b ■ V y (b ■ V y u)) ■ VyV dy. 

JR m JR m 

□ 

The matrix fields F G dom(L 2 ) and E = L(F) G dom(L) have the following properties. 

Proposition 5.1 For any u,v G C 1 (M m ) which are constant along the flow of b we have in 
V'(R m ) 

DVyU ■ VyV — {D)q VyU • V yV = — 6 ■ V y{EV yU ■ V yV) = — diV ?/ (6 <g) bVy(FVyU • V yV)) 
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and 

{EVyU ■ VyV) = (FVyU ■ V ' yV) = 0. 

In particular 

/ EV y u ■ V y v dy = (EV y u ■ V y v) dy = 
jR m Jm. m 

/ FVyU ■ V y v dy= (FVyU ■ V y v) dy = 
saying that (dw y (EV y u)) = (div y (FV y u)) = in V '(R m ). 

Proof. Consider ip G C\ (R m ), u, v G C 1 (M m ) such that u s = u, v s = v, s G R and the matrix 
field U = (pQ~ 1 V y v <g> V y uQ~ l G Hq. Actually U G dom(L) and, as in the proof of the last 
statement in Proposition I3.13[ one gets 

L(U) = (b ■ V y ip)Q~ l V y v <g> VyuQ" 1 + <p L(Q~ l V y v <8> VyiiQ- 1 ) 
= (b ■ Vyif)Q~ 1 V y v ® VyuQ' 1 

since Q~ x V y v (g> VyuQ^ 1 G ker(L). Multiplying by E7 the equality D - (D) Q = -L(E), 
E = L(F), one gets 

/ <p(D-(D) Q )V y u-V y vdy = -(L(E),U) Q = (E,L(U)) Q = ( (b ■ V y f>){EV y u ■ V y v) dy 

implying that DV y u ■ V y v = (D) Q V y u ■ V y v - b ■ V y (EV y u ■ V y v) in V '(R m ). Multiplying 
by U the equality E = L(F) yields 

/ <pEV y u-Vyvdy = (E,U)Q = (L(F),U)Q = -(F,L(U)) Q = - [ (&• V y (p)FV v u- V y v dy. 
We obtain 

EVyU ■ VyV = b ■ Vy{FVyU ■ V ' yV) UL V '(R™) 

and thus 

DVyU ■ VyV - (D) Q VyU • VyV = ~6 ■ Vy(EVyU • V>) = ~6 ■ Vj / (6 ■ Vy{FVyU • Vj,«)) 

in P'(R m ). Consider now J7 = </<) _1 V y 77 <g> VyuQ^ 1 with 99 G ker(6 • V y ). We know that 
L(C7) = and since, by construction F G (kerL)- 1 , we deduce 

/ LpFV y u ■ V y v dy = (F, U)q = 

saying that {FV y u ■ V y v) = 0. Similarly E = L(F) G (ker L) 1 - and {EV y u ■ V y v) =0. a 
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Remark 5.1 Assume that there is Uq (eventually multi-valued) satisfying Uo(Y(s;y)) = 
u o(y) + s > s € M.,y G M m . Its gradient changes along the flow of b exactly as the gra- 
dient of any function which is constant along this flow cf. Remark \3. 1\ We deduce that 
Q v \l y v <%> VyuQ" 1 G kerL for any u, v G ker(6 • V y ) U {uq} and therefore the arguments in 
the proof of Proposition I5.il still apply when u, v G ker(6 ■ Vy) U {uo}- I n the case when m — 1 
independent prime integrals {ui, u m —\] of h are known, the matrix fields E,F come, by 
imposing for any i,j G {0, 1, m — 1} 

— b ■ Vy(EV y Ui ■ VyUj) = DVyUi ■ VyUj — (DVyUi ■ VyUj) , (EVyUi ■ y y Uj) = 

and 

b ■ V y (FV y Ui ■ VyUj) = EVyUi ■ VyUj, (FV y Ui ■ V y Uj) = 0. 
We indicate now sufficient conditions which guarantee that the range of L is closed. 

Proposition 5.2 Assume that (|13p . (|14p , (|28p hold true and that there is a matrix field R(y) 
such that ()24p holds true. Then the range of L is closed. 

Proof. Observe that (|24p implies (|17p . Indeed, it is easily seen that b ■ V y R + Rd y b = in 
V'(R m ) is equivalent to R = R s d y Y(s; •), s G M. We deduce that P = iT 1 t RT l satisfies 

G{s)P = d y Y-\s; .)P B %Y-\a\ ■) = d y Y-\s; -)RJ X 'R' 1 ^yY'^s; •) = R' 1 'iT 1 = P 

saying that [b, P] = in T>'(R m ). Therefore we can define L as before, on Hq, which 
coincides in this case with {^4 : RA t R G L 2 (R m )}. We claim that i o L = (b ■ W y ) o i where 
i : Hq — > L 2 (R m ), i(A) = RA l R, A G Hq, which comes immediately from the equalities 

(i o G(s))A = RG(s)A t R = Rd y Y- x (s; ^AjdyY'^R = i? s A s *i? s = (i(A)) s , s € R,i G Hq. 

In particular we have 

kerL = {A G Hq : i(A) G ker(6 • V y )} 

and 

(kerL) -1 = {^4 G Hq : [ i(A) : U dy = V U G ker(6 • V y )} 

= {AGH Q : i(A) G (ker(&- Vy)) 1 }. 

For any >1 G (kerL)- 1 we can apply the Poincare inequality (|28j) to «(j4) G (ker(6 • V^)) and 
we obtain 

\A\ Q = \i(A)\ L * < C P \b ■ V y (i(A))\ L 2 = C P \i{L{A))\ L , = C p \L{A)\q. 
Therefore L satisfies a Poincare inequality as well, and thus the range of L is closed. 
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Remark 5.2 The hypothesis b ■ W y R + Rd y b = in £>'(IR m ) says that the columns of R 1 
form a family of m independent vector fields in involution with respect to b, cf. Proposition 



R- 1 (y) = d y Y(s-y)R- l (y) 1 s£l, y £ R m . 

Remark 5.3 For any U £ kerL, that is i(U) € ker(6 • Vy) ; we have 

[ R(D - (D) Q fR : i(U) dy = 0. 

As (D)q € kerL, we know that i((D)q) = R (D)q t R € ker(fo • Vy) and thus the matrix field 
R{D}q t R is the average (along the flow ofb) of the matrix field RD t R, which allows us to 
express (D)q in terms of R and D 

R(D) Q t R= (RD t R). 

From now on we assume that (|24p holds true. Applying the decomposition of Theorem 12.31 
with the dominant term u € ker(6 • V y ) in the expansion and any v £ Cf (R m ) yields 

f (D - (D) Q )VyU -VyVdy=- f FVyU ■ Vy(b ■ Vy{b • V yV)) (fy. 

From ([58]) one gets 

f (D - (D) Q )V y u ■ V y v dy - [ u l b-Vy(b-V y v)dy = 

it™ Jm. m 

and thus 

u 1 = diVy(FVyu) + V 1 , v 1 Gker(6- V y (b- Vy)) =ker(6- V y ). (59) 

Notice that (n 1 ) = v , since {div y (FV y u)} = 0, cf. Proposition 15.11 The time evolution for 
v 1 = (u 1 ^ comes by averaging (|57|) 

dtv 1 - (divyiDVyV 1 )) - {dwy(DVy{div y (FV y u)))) = 0. 

As v 1 £ ker(6 • Vy) we have 

- (diVyODVX )) = -diV y ((D) Q VyV 1 ) 

and we can write, with the notation w 1 = div y (F\7 y u) 

d t {u + eu 1 } - div y ({D) Q V y {u + en 1 }) = ed t w l - edw y ((D) Q VyW 1 ) + e (divyiDVyW 1 )) . 

(60) 
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But the time derivative of w is given by 

dtw 1 = dw y (FVyd t u) = dW y (FV y (div y ({D) Q V y u))) 

which implies 

dtw 1 - div y ((D) Q Vyw 1 ) = div y (FV y (div y ({D) Q X7 y u))) - div y ({D) Q V„(div v (FV v «))) 

= -[6iv y ((D) Q V v ),diVy(FV y )]u. 

Up to a second order term, the equation (|6U|) writes 

d t {u + eu 1 } - dw y ((D) Q V y {u + eu 1 }) + e[div y ({D) Q V„), div y (FV y )]{u + eu 1 } 

- e (div y (DV y (div y (FV yU )))) = 0(e 2 ). (61) 

We claim that for any u € ker(6 • V y ) we have 

(dw y (DV y (dw y (FV y u)))} = {div y (EVy(divy(EVyu)))) . (62) 

By Proposition 15.11 we know that (div J/ (i ? Vj / u)) = 0. As L({D)q) = we have 

[b ■ V y , -div y ((D) Q Vy)) = -div y (L((D) Q )Vy) = 

and thus div y ((D)q Vy) leaves invariant the subspace of functions which are constant along 
the flow of b. By the symmetry of the operator divy ((.D)q S/ y ), we deduce that the subspace of 
zero average functions is also left invariant by div y ({D)q Vy). Therefore /div y ((Z))g V y (div y (FV 
and 

(dw y (DV y (dw y (FV y u)))) = (div y ((D - (D) Q )Vy(div y (FVyu)))) . 
Thanks to Theorem 12.31 we have 

div„((D - (D) Q )V y ) = [b ■ V y , [b ■ Vy, -dw y (FV y )} } 

= [b-Vy,-dhy(L(F)Vy)} 
= [b-Vy,-diVy(EVy)) 

which implies that 

(dwy(DVy(dw y (FVyu)))) = (div y ((D - (D) Q ) Vy (div y (FV y u) ) ) ) 

= {diV y (EV y (b ■ V y (diVy(FV y U)))) - b ■ Vy(diVy(EVy(dWy(FVy U ))))) 
= (diVy(EVy(b ■ V y (dw y (FVyU))))) . 
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Finally notice that 

-divy(EVyu) = -divy(L(F)V y u) = [b ■ V y , -div y (FV y u)] = -b ■ V y (dw y (FS/ y u)) 

and ([62]) follows. We need to average the differential operator div y (E'Vy(divy(E\/ y ))) on 
functions u E ker(6 • V y ). For simplicity we perform these computations at a formal level, 
assuming that all fields are smooth enough. The idea is to express the above differential 
operator in terms of the derivations t R which commute with the average operator (see 
Proposition 13. 6|) . since the columns of R^ 1 contain vector fields in involution with b(y). 

Lemma 5.1 Under the hypothesis (|24|) . for any smooth function u{y) and matrix field E(y) 
we have 

divy(EVyu) = div y (R t E) • (*i? -1 V y «) + RE t R : (*fl -1 V tf <8 'iTV^K (63) 

Proof. Applying the formula divy(A^) = div^-A • £ + * A : dy£, where A(y) is a matrix field 
and £(y) is a vector field, one gets 

divy(EV y u) = div y (E l R t R- x V y u) = d\v y {R t E) ■ (*i? _1 V y «) + R t E: S^(*fl -1 V y u). 

The last term in the above formula writes 

R l E : d y ^RT l V y u) =R t E t R 'iT* : S^(*fl -1 V v u) 
= R t E t R: dy^R~\yU)R- 1 
= RE t R: l R~ l tdy^R^Vyu) 
= RE t R: CR^Vy ® t R~ 1 V y )u 

and ([63]) follows. | — | 

Next we claim that the term (drvy(E'Vy(diVy(ffVyu)))) reduces to a differential operator, if 
u G ker(6 • V y ). 

Proposition 5.3 Under the hypothesis (|24p . /or any smooth matrix field E there is a linear 
differential operator S(u) of order four, such that, for any smooth u € ker(6 • V y ) 

(div v (^„(div tf (SV tf u)))) = S(u). (64) 
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Proof. For any smooth functions u,<p G ker(6 • V^) we have, cf. Lemma 15, II 

/ (divy (EV y (divy (EV y u)) ) } ip dy = / div^ (E V y (div y (EV y u) ) )tp dy 

JR m JR m 

= / d\v y {EV y u) d\Vy{EV ytp) dy 

JR m 

= [ {div y {R t E) ■ CR^Vyii) + RE t R: ( t iJ -1 V v <8> *R~\y)u} 
Jr™ 

x {dw y (R l E) ■ fR^Vyip) + RE t R: (*i? _1 V y ® *R -1 V„M dy 
= / [div„(12 *S) ®div„(12 : V y u ® V„p] dy 

+ / [RE *i? <g> div„(J2 % E)] : [C-R -1 V tf ® *i? -1 V v )« <8> t fl" 1 V tf ^] dy 
Jr™ 

+ / [div„(/i iLE *E] : [C-R^Vyu) ® (*fJ -1 V v *22 -1 V tf )¥?] dy 

+ / [EE tR®RE t R]:[( t R- 1 V y ® t R~\ y )u®( t R- 1 V y ® t R~ 1 V y )<p]dy 

jR m 

Recall that t R~ 1 S/ y leaves invariant ker(6 • Vy) and therefore 

t R~ 1 VyU <g) t R~\ y (p G ker(6 • V y ) 

implying that 

/ [div a (i? *E)®div tf (i2 *£?)] : ^BT 1 V y u ® t ET 1 V v tp] dy 
= / (div y ( J R *E) Odivj / ( J R *£)} : ^R- x V v u® t Er\ y <p] dy. 

jR m 

Similar transformations apply to the other three integrals above, and finally one gets 
/ {divy(EVy(dw y (EVyu)))} ip dy = / X : [V R u <g> V R ip] dy 

JR m JR m 

+ Y : [{V R ® V R )u <g> V'V] dy 
+ Z: [V R u <g> (V R <g> V fl )v?] dy 
+ I T : {(V R ® V K )u ® (V fl ® V K )v9] dy 

JR m 

= h(u,tp) +I 2 (u,(p) +I 3 (u,(p) +h(u,<p) 

where V R := *.R~ V„ and X,Y,Z,T are tensors of order two, three, three and four respec- 
tively 

X tj = (dw y {R *E)i div y (R t E) j ) , i,j G {1, ...,m} 
y ijfc = ((RE */J) fj div y (i? , Z yfc = (div y (# (-RE *i?) ifc ) , i, j, k G {l,...,m} 

T ijH = ((RE i, j, A;, / G {1, m}. 
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Integrating by parts one gets 

I^u, ^) = f XV R u ■ V R ip dy= [ R- 1 XV R u ■ V y ip dy = f Si(u)tp dy 
Jm m Jm m jR m 

where S\(u) = —div y (R~ 1 X'V R u). Notice that the differential operator 

£ -> divyiR-^) = div y { 'iT 1 ) ■ £ + : dyi 

maps (ker(6 • V y )) m to ker(6 • Vy), since the columns of R~ l contain fields in involution with 
b, and therefore S% leaves invariant ker(6 • Vy), that is, for any u E ker(6 ■ V y ), £ = XV R u E 
(ker(6 • V y )) m and Si{u) = -div^iT^V^u) = -div^iT 1 ^ E ker(6 • V„). Similarly we 
obtain 

h(u,ip) = / S 2 (u)<pdy, I 3 (u,(p)= / S 3 (u)(pdy, h(u,<p) = / S^u)^ dy 

JR m JR m J]R m 

where 82,83,84 are differential operators of order three, three and four respectively, which 
leave invariant ker(6 • S/ y ). We deduce that 



(dwy(EX7y(div y (EV y u)))) ip dy = / S(u)(p dy 
for any u, (p E ker(6 • V y ), with S = S\ + S2 + S3 + S4, saying that 

(div y (EVy(divy(EVyu)))} - S{u) _L ker(6 • V„). 
But we also know that 

(divy(E\7 y (dwy(EVyu)))} - S(u) E ker(6 • V y ) 

and thus (fM|) holds true. 
Combining (IBB . (IS2 L we obtain 



□ 



^{u + en 1 } - div w ((D> g V y {u + en 1 }) + e[div y {{D) Q V y ), div y (FS/y)]{u + eu 1 } 

-eS{u + eu l ) = 0{e 2 ) 

which justifies the equation introduced in (|25p . The initial condition comes formally by 
averaging the Ansatz Q 

{u £ ) = u + ev 1 + 0(e 2 ). 



One gets 



v (0, •) = w- lim 



e\o e 

implying that n 1 (0, •) = Vi n + div^-FVyUm), cf. (f59|) . which justifies fl26j) , 
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6 An example 

Let us consider the vector field b(y) = y := (y2,—yi), for any y = (yi , j/2) £ R 2 and the 
matrix field 

( A i(y) \ 

D(y)= , ytR 2 

\ X 2 (y) ) 

where Ai,A2 are given functions, satisfying min ygK 2{Ai(y), A2(y)} > d > 0. We intend to 
determine the first order approximation, when e \ 0, for the heat equation 

d t u £ - divy(D(y)V y u £ ) - ^div y (b(y) (g> b{y)V y u £ ) = 0, (t, y) G R + x M 2 (65) 

with the initial condition 

u £ {o,y) = Uin {y), yei 2 . 

The flow of b is given by y) = 1Z(—s)y, s£l,j;£M 2 where 7£(a) stands for the rotation 
of angle odt The functions in ker(6 • V^) are those depending only on \y\. Notice that the 
matrix field 



\y\ 



yi V2 



satisfies b ■ V y R + Rd y b = and Q = t RR = I 2 . The averaged matrix field {D)q comes, 
thanks to Remark |5. 3\ by the formula R (D)q l R = i^RD and thus 



(D) Q = t R(RD t R)R, {RD l R) 



Ai3/|+A 2 yf \ / (Ai-A 2 )j/ii/2 

\yp / \ \vr 

(Ai-A 2 )yii/2 \ / Aiy^+A 2 y| 
\yp / \ W 



In the case when Ai, A2 are left invariant by the flow of b, that is Ai, A2 depend only on \y\ 
it is easily seen that 

(w) = (w) = r (w) = 

and thus 



x R 2 



(D) Q =tR±±^I 2 R=±±^I 2 . 
The first order approximation of (|65p is given by 

dtU _ diVj/ ( M+*aM Vy tt) = 0, (t, y) e 
u(0,y) =u in (y), y G R 2 . 

We consider the multi-valued function uo(y) = —9(y), where y = |y|(cos 9(y), sin 6(y)), which 
satisfies b ■ V y uo = 1, or uo(Y(s; y)) = uo(y) + s. Notice that the averaged matrix field {D)q 
satisfies (with u\{y) = \y\ 2 /2 G ker(6 • V y ) ) 

VyUi ■ (D)q VyUj = (VyUi • DVyUj) , l,j G {0, l} 

as predicted by Remark 13.71 
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A Proofs of Propositions [33} [33], ULS, HL2 



Proof, (of Proposition 13, 4|) For simplicity we assume that b is divergence free. The general 
case follows similarly. Let c(y) be a vector field satisfying (|29p . For any vector field <fi G 
(R m ) we have, with the notation u T = u(Y(t; •)) 



c • -<p) dy= (c h - c) ■ <f) dy = / (d y Y(h; y) -I)c-(j> dy. 

n JR m JR m 

Multiplying by /i -1 and passing to the limit when h — > imply 

c(6 • V y </>) dy = dybc ■ 4> dy 
n JR m 

and therefore (6 • V y )c - <9 y 6c = in V '(R m ). 

Conversely, assume that [b, c] = in V^W 71 ). We introduce e(s,y) = c(Y(s;y)) — 
d y Y(s;y)c(y). Notice that e(s, •) G L^IR" 1 )^ G R and e(0, •) = 0. For any vector field 
<\> G Cl(R m ) we have 



E^s) := / e(s, y) • <j){y) dy = c(y) ■ dy - d y Y(s; y)c{y) ■ <j)(y) dy 

JR m JR m JR m 

and thus 

^-E^s) = - [ c(y) ■ ((b ■ V„)0)_ a dy - / d,(6(Y( S ; y))) c(y) • <f>(y) dy 

as JR m JR m 

= - [ c ■ (b ■ V y )^ s dy - [ d y b(Y(s;y))d y Y( S ;y)c(y)-<P(y)dy 

JR m JR m 

d y b c(y) ■ 4>- s dy- I d y b(Y(s; y))d y Y(s; y)c(y) ■ 4>{y) dy 

n JR™ 

d y b(Y(s; y))(c(Y(s; y)) - d y Y(s; y)c(y)) ■ </>(y) dy 

n 

e(s,y)- t d y b(Y( S ;y)) ( t } (y) dy. 

n 

In the previous computation we have used the fact that the derivation and tranlation along 
b commute 

((&-V„)0)_ a = (&• V„)#_ 5 . 
After integration with respect to s one gets 

M s ) = I* f V) ■ %b(Y(r; y))<f>(y) dy dr. 
JO JR m 

Clearly, the above equality still holds true for any <f> G C c (M. m ). Consider R > 0,T > and 
let K = || dyb o Y\\ L ooi\_rp,T]xB R )- Therefore, for any s G [— T,T] we obtain 

\\e(s,-)\\ L ^(B R ) = sup{| J B 9i (s)| : (f> G C c (Br), |H| L i (R m) < 1} 



< K 



|e(r, •)IU-(B J ,)dr 
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By Gronwall lemma we deduce that ||e(s, Olli/WS^) = for — T < s < T saying that 
c(Y(s; y)) - d y Y(s; y)c(y) =0,s6l,i/G M m . □ 



Proof, (of Proposition 13, 5p 

1. =^ 2. By Proposition 13.41 we deduce that c(Y(s;y)) = d y Y(s;y)c(y) and therefore 
/ (c • Vyu)v. s dy= f c(Y(s; y)) ■ (V y u)(Y(s; y))v{y) dy 

JR m JR™ 

c(y) ■ 'OyYis; y)(V y u)(Y(s; y))v(y) dy = [ (c(y) ■ V y u s )v(y) dy. 



2. ==> 3. Taking the derivative with respect to s of (|30p at s = 0, we obtain (|3ip . 3. ==> 1. 
Applying ([3T]) with v G C^(R m ) and m = ynp{y), f G C%(R m ), <p = 1 on the support of v, 
yields 

/ a b ■ V y v dy + c • Vyh v(y) dy = 



saying that b ■ V y Ci = (d y b c)i in V '(R m ), i G {1, m} and thus [6, c] = 6 • V y c — c^frc = 

□ 



in V(R m 



Proof, (of Proposition I3.8|) The arguments are very similar to those in the proof of Proposi- 
tion 13.41 Let us give the main lines. We assume that b is divergence free, for simplicity. Let 
A(y) be a matrix field satisfying (|33p . For any matrix field U G C^(M m ) we have 

/ A(y):(U(Y(-h;y))-U(y))dy= [ (A(Y (h; y)) - A(y)) : U (y) dy 

JR m JR m 

r (d y Y(h;y)A(y) t d y Y(h;y) - A(y)) : U(y) dy 
{(d y Y(h; y) - I)A(y) t d y Y(h; y) : U(y) + A{y)\d y Y{h- y) - I) : U(y)} dy. 
Multiplying by 4 and passing /i^Owe obtain 

' A(y) : (b ■ V y U) dy= f (d y bA(y) + Aiyfdyb) : U(y) dy 



saying that [6, A] = in P'(R m ). 

For the converse implication define, as before 

f(s,y) = A(Y(s;y))-d y Y(s;y)A(y) t d y Y(s;y), s£l, y G M. m . 

For any G C^(R m ) we have 

%(«):= / f(s,y):U(y)dy 

m 

A(y) : C/(Y(- S ; y)) dy - [ d y Y(s; y)A(y) t d y Y(s; y) : U(y) dy 

n JR m 
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and thus 

^-Fu(s) = - [ A(y) : ( (b ■ V y )U )_ a dy - f d y (b(Y(s; y)))A{y)%Y{ a] y) : U(y) dy 

as JR m JR m 

d y Y(s;y)A(y) t d y (b(Y(s;y))) : U(y) dy 

n 

A(y) : (b ■ V y )U- s dy - f d y b(Y(s; y))d y Y(s; y)A(y) t d y Y(s; y) : U dy 

n JR m 



d y Y(s;y)A(y) t d y Y(s;y) t d y b(Y(s;y)) : U(y) dy 

n 

{d y b(Y(s; y))f(s, y) + f(s, y)*^6(y(a; y))} : U(y) dy 

■n 

= [ f(s, y) : {%b(Y(s; y))U(y) + U(y)d y b(Y(s; y))} dy. 

The previous equality still holds true for U € C c (M m ), and our conclusion follows as in the 
proof of Proposition 13.41 by Gronwall lemma. 

Proof, (of Proposition I3.9() 

1. =^ 2. By Proposition ES] we deduce that A(Y(s;y)) = d y Y(s;y)A(y) t d y Y(s;y). Using 
the change of variable y — > Y(s; y) one gets 

A{y)V y u-Vyvdy= [ A(Y( S ;y))(V y u)(Y(s;y)) ■ (V y v)(Y(s;y)) dy 

n 

A( y )%Y(s;y)(y yU )(Y(s;y)) • t Wa;y)(V w u)(y(s;y)) dy 

n 

A(y)V y u s ■ V y v s dy. 

n 

2. ==>■ 3. Taking the derivative with respect to s at s = of the constant function s — > 
f Rm A(y)V y u s ■ V y v s dy yields 

/ A(y)V y (b- V y u) ■ V y vdy+ I A(y)V y u-\7 y (b-Vyv)dy = 0. 

jR m JR™ 

3. ==>■ 2. For any u, v G C^(M m ) we can write, thanks to 3. applied with the functions u s ,v s 

-y- / A(y)V y u s ■ V y v s dy = f A{y)V y { (b ■ V y u) s ) ■ V y v s dy 

& s JR m JR m 

+ / A(y)V y u s ■ V y ( (b ■ V y v) s ) dy 

JR m 

= / A(y)V y (b ■ V y u s ) ■ V y v s dy 

JR m 

+ / A(y)V y u s ■ V y (b ■ V y v s ) dy = 0. 

JR m 

Therefore the function s — > j Rm A(y)V y u s ■ V y v s dy is constant on E and thus 
/ A(y)V y u s ■ V y v s dy = A(y)V y u ■ V y v dy, seR. 

jR m jR m 
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Up to now, the symmetry of the matrix A{y) did not play any role. We only need it for the 
implication 2. ==> 1. 
2. => 1. We have 

/ A(y)V y u ■ V y v dy = A(y)V y u s ■ V y v s dy 

A{y) t d y Y{s; y)(V y u) s ■ y)(V y v) s dy 

d y Y(s;y)A{y) t d y Y(s;y)(V y u) s • (X7 y v) s dy 

n 

= [ {d y YA t d y Y)_ s V y u-V y vdy 

where {d y YA t d y Y)^ s = d y Y(s;Y(-s;y))A(Y(-s;y)) t d y Y(s;Y(-s;y)). We deduce that 

/ (A(y) — (d y YA t d y Y)^ s )'V y u ■ V y v dy = 0, u,v£Cl(R m ). 

Since A{y) — (d y YA t d y Y)^ s is symmetric, it is easily seen, cf. Lemma IA.1I below, that 
A{y) - {d y YA tdyY)^ = 0. Therefore we have A(Y(s;y)) = dyY(s;y)A(y) t dyY(s;y) : s G 
R, y G M. m and by Proposition ES] we deduce that [6, A] = in V '(R m ). n 

Lemma A.l Consider a field A(y) G Lj^W 71 ) of symmetric matrix satisfying 

[ A(y)V y u ■ V y v dy = 0, u, v G Cl{R m ). (66) 
Therefore A{y) = a. a. y G M. m . 

Proof. Applying ([66]) with Vj = yjv, v G (M m ), U{ = yiip(y) where <p G C^(M. m ) and <p = 1 
on the support of v, yields 

/ A(y)e i -(y j V y v + ve j )dy = 0. (67) 

JM. m 

Applying ([66]) with v and Uij = yiyjip(y) one gets 

/ A(y){y j e i + y i e j )-V y vdy = 0. (68) 
Combining (|67|) . (|68p we obtin for any i,j G {1, ...,m} 

2 / (A(y)ei ■ ej) v(y) dy = (A(y)e< • + A(y) ej • ei)u(y) dy = 
saying that A(y) = 0, a.a. y G M m . q 
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